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ABSTRACT 

We have identified an optical/X-ray binary with orbital period Pb = 5.47h as the likely counterpart of the 
Fermi source 2FGL J2039.6—5620. GROND, SOAR and DES observations provide an accurate orbital period 
and allow us to compare with the light curve of an archival XMM exposure. Like many short-period optical 
X-ray binaries associated with FAT sources this may be a interacting (black widow/redback) millisecond pulsar 
binary. The X-ray light curve is consistent with the emission associated with an intrabinary shock. The optical 
light curve shows evidence of companion heating, but has a peculiar asymmetric double pei. The nature of this 
optical structure is not yet clear; additional optical studies and, especially, detection of an orbital modulation in 
a 7 -ray pulsar are needed to elucidate the nature of this peculiar source. 

Subject headings: gamma rays: stars — pulsars: general 


1. INTRODUCTION 

As identification of the 1873 O.l-lOOGeV sources in 
the second Fermi Large Area Telescope (FAT) catalog 
(iNolan et al.l |2m^ progresses, blazars and s pin-powered 
pulsa rs increasingly dominate the associations (lAcero et al.l 
l2015h . For the brightest remaining unidentified sources 
gamma-ray spectral curvature and variability provide an 
excellent disc rimination between these two possibilities 
(lRomanill2012h . 2FGL J2039.6-5620 = 3FGL 2039.6-5618 
(hereafter J2039) is a relatively bri ght fn.i -1 nn OrV = 1-7 ± 
0.1 X lO”'^'^ ergcm“^ (25.4g, lAcero et all 120151) source 
with a good localization at Galactic latitude -37°, a 3FGL 
“variability index” value of 34, and a “curvature significance” 
of 5.1 ct. Thus it is one of the least variable bright unidenti¬ 
fied sources and displays significant spectral curvature, plac¬ 
ing it well within the pulsar zone of the curvature-variability 
plan e. Accordingly it h as been repeatedl y searched for ra¬ 
dio (ICam ilo et al. Tl2015h and gamma-ray (lAbdo et ^ 120091: 
iPletsch & Clarkll2014ft pulsations. The lack of any detection 
is best understood if the source is a millisecond pulsar rapidly 
accel erated and/or wind-shrouded in a close, interacting bi¬ 
nary d^ mani & ShawmOlTt !Romanil[2012l: Ide Martino et al.l 
l2^lRov et al.ll2015h. 

We have examined new and archival optical and X-ray data, 
discovering a short period (5.47 h) binary. The system shows 
evidence of heating, and the X-ray light curve is reminis¬ 
cent of other shrouded (black widow BW and redback RB) 
millisecond pulsar binaries. The optical light curve is oddly 
double-peaked and asymmetric. Nevertheless we conclude 
that this is another short period gamma-ray pulsar binary. 
With our well measured period and position, this is a fruit¬ 
ful target for renewed pulse searches. 

2. PHOTOMETRY AND ORBITAL-PERIOD ESTIMATE 

J2039 has a good 3FGL localization (2.7' radius) and we 
initially observed this field as part of our FAT Unidenti¬ 
fied source campaign, using the Goodman High Through¬ 
put Spectrograph (GHTS) at the 4.2 m SOAR telescope on 
2013 August 10 (MJD 56514, UT dates are used through¬ 
out). A sequence of 5 x {g'r'i') 180 s images plus a 600 s 
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Fig. I.— Stacked SOAR/GHTS g' image of the J2039-5617 field. 

Ha frame were obtained, spaced over 2 hours and largely 
covering the error ellipse. Despite relatively poor condi¬ 
tions, we were able to search for variability for counter¬ 
parts as faint as g ~ 22. One star (Figure 1), coincident 
with an X-ray source in archival data (initially SWIFT, later 
XMM exposure, see below) was found to vary in a corre¬ 
lated manner in all three bands, suggesting a characteristic 
period ^ 3 — 6 h. The USNO-BFO position of this source is 
a = 20^39“34.987", (5 = -56°17'09.04" (2000.0, aRA = 
0.06", (Jdec = 0.13"). The catalog gives a proper motion 
estimate /ira = 14±4masy“'^ and /tra = —16±9masy“'^. 
No Ha structure was found in the field. 

We searched for archival observations in this direction 
and f ound an excellent series of GROND (iGrenier et alJ 
l2008h exposures of the J2039 error ellipse, with 72 x 115 s 
p'r'i'z'(H-JHK) exposures each on 3 days Aug 16-18 2014 
(MJD 56885-56887). Measuring these optical frames, we ob¬ 
tain the light curves shown in figure 2. This data set showed 
a clear double-peaked light curve, with asymmetric maxima 
and minima and an orbital period Pr = 0.22815 ± 0.00015d 
(5.47 h). The colors vary slightly with evidence for heating 
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GROND error bars plotted on the first cycle only. GHTS (circle) and DES 
(triangle points are also shown. The DES z‘ offset may be due to the sub¬ 
stantial differences in this filter. Phase markers on the second period of the 
g' curve emphasize that the second maximum and minimum are eai'ly by 
0.05. A single period of the g' — r' color is also shown, on a 
magnified scale. 

(minimum g' — r') at optical maximum. 

The quest for direct detection of blind pulsations is greatly 
helped by a precise orb ital period, which reduces the ac¬ 
celeration search space (iPletsch et al.l I20l3) . Accordingly, 
we compared with light curve phasing of our GHTS obser¬ 
vations, 371 d earlier. This refined the period to 0.228116 
±V(1.4 X 10“"^) d where the aliases about the best minimum 
are set by orbit slips between the two observation sets. 

The target was covered in seven exposures in the DES early 
data release, in two frames each of g', r' and i', and a single z' 
frame. We performed relative photometry against a set of the 
field stars from the GROND images, to put the magnitudes 
on a consistent scale. The observation dates ranged from 
MJD 56538 to MJD 56558, between our other two epochs. 
Accordingly these points break the alias degeneracy, clearly 
preferring N = 0, the primary minimum. We thus estab¬ 
lish Pb — 0.228116 ± 0.000002, with barycentric epoch for 
optical maximum MJD 56885.085. As discussed below it is 
unclear how this photometric epoch relates to the orbital kine¬ 
matic epoch, i.e. Tasc- 

In Figure 2 the GHTS and DES points are seen, calibrated 
against the GROND flux scale, showing good agreement with 
the mean light curve. We also plot the DES z' point, which 
lies —O.lSmag from the GROND z' curve. This discrepancy 
may be cross-calibration error (GROND z' is split between 
DES z' and Y). The DES data are however very spread in 
time; the z' point was 2.5 orbits off any other observation. 
The second g' point is also ^ 3a from the expected flux 
and was obtained 5 d before the first DES g' point. If, like 
many other BW/RB systems, the companion undergoes occa¬ 
sional flares, this can explain these discrepancies. However 
the small spread in GROND photometry over 3 d and the oth¬ 
erwise good match of the GHTS and DES points imply that 
the source is generally stable to better than ~ 0.02mag. 

J2039 is also covered in the Catalina Sky Survey (CSS; 
http://www.lpl.arizona.edu/css/index.htmli Drake et al. 2009) 


photometry archive with data extending 9.0 y before the 
GROND exposures. These unflltered magnitudes with a typ¬ 
ical error ^ O.Smag are too shallow to usefully probe the 
light curve modulation. However, they do serve to confirm 
the source’s general quiescence; out of 217 observations, only 
one point with am < 0.4mag had m < 18 (likely a measure¬ 
ment error in any event). 

The most remarkable aspect is of this light curve are the two 
asymmetric maxima and minima. The first minimum follows 
the large maximum by ~ 0.25, but the second maxi¬ 
mum and minimum arrive « 0.05 too early to lie on the 
opposite side of the orbit (dotted lines in Figure 2). The opti¬ 
cal color is bluest at optical maximum, but plateaus across the 
first minimum and second maximum, decreasing to its reddest 
value at the second minimum. Clear there is off axis structure 
in the system, with heating at two (nearly) opposite poles or, 
possibly an occultation at the phase of the first minimum. 



bined MOS1+MOS2. = 0 is at optical maximum as in Figure 2. 

3. ARCHIVAL X-RAY LIGHT CURVE 

Archival SWIFT exposures show a few counts from this 
counterpart, but the X-ray source is best measured in a 41.6 ks 
live time XMM-Newton exposure (ObsID 0720650301, PI: 
Mignani) on MJD 56575. This data set clearly shows the op¬ 
tical counterpart as the brightest source in the 3FGL error el¬ 
lipse. With an observation between the GROND and GHTS 
epochs, the predicted orbital phase is known to S(j)B = 0.012 
dominated by the period uncertainty. Folding the X-rays, we 
generated the PN and combined MOSl/2 light curves (figure 
3). The peak is coincident with the optical maximum. The 
double structure, while of modest statistical significance, is 
common to both instruments. 

The source is unresolved to XMM and delivered 0.024 cps 
in 37 ks effective PN camera exposure. A power law spec¬ 
tral fit determines a photon index E = l.Sd^'^Q'Jg and an up¬ 
per limit on the absorption of Nb < 1.9 X lO^Ocm-^ (90% 
errors). The corresponding unabsorbed 0.5-10keV flux is 
8.2lg'5 X 10“^^ergcm“^s“^. If a thermal component (e.g. 
from an underlying neutron star) is present a larger Nh might 
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be allowed. However present data do not require such a com¬ 
ponent and the absorbed power law fit is quite adequate, with 
x2/DoF=1.15. 

4. SYSTEM MODELING AND CONCLUSIONS 

The source bears many similarities to several other short- 
period shrouded pulsar binaries. At Pb — 5.47h it may be 
either a redback or a black widow. However, at this period 
it is much closer to the shrouding line in the former case and 
so it is more likely a redback, given the lack of detection in 
repeated radio searches. Although the system properties will 
remain quite uncertain until a pulse detection is made, we can 
draw some tentative conclusions from the binary data. 

Our X-ray spectral fit indicates very little extinction, with 
the Nh up per limit correspond ing to Ay = 0.07 for the con¬ 
version of iFoight et al.l (l2015h. This is about half the tota l 
Ay = 0.16 in this direction (ISchlaflv & Finkbeinedl20ll]) . 
At this orbital period a main sequence companion would fill 
its Roche lobe for Me = 0.61Mq so we take this as an up¬ 
per limit to the companion mass. Such a star would have 
Te// = 4180K and g' — r' = 1.27, so the observed color 
at minimum g' — r' = 0.78 indicates appreciable heating. In 
fact, the colors are even bluer at maximum with g' — r' = 0.73 
equivalent to Te// = 4980 K. The origin of these increased 
temperatures is unclear; the companion may be partly degen¬ 
erate or it may experience tidal heating. However in analogy 
with other shrouded MSP, it is likely that the effect of pulsar 
heating, possibly mediated by an intrabinary shock, is domi¬ 
nant. If we consider just the temperature rise to maximum, we 
derive an irradiation temperature ~ 4200 K and for an orbital 
separation a = 2.O(Mtot/2M0)^/^i?0, an isotropic irradia¬ 
tion luminosity 4.2 X 10^^ergs“^. Lower companion masses 
require additional heating power, although the relatively large 
TMin implies a low inclination angle or, possibly, rapid trans¬ 
port of heating flux to the companion night side, a s invoked 
for some hot Jupiters (e.g. [Heng & Showmanl2015h . 

This heating flux is a lower limit to the true pulsar output. 
iBreton et al.l (120 13h infer a typical heating efficiency p = 0.1 
(some BW/RB appear to have larger heating efficiencies), 
which would imply a spin-down power E k, 4xl0^^ergs 
Another power estimate can b e derived from the gamma-ray 
flux 1.7 X 10~^^erg cm~^s~^: lAbdo et al.l (l2010h have noted 
a heuristic gamma-ray luminosity T-y = (10^^ergs“^T)^/^, 

which gives us a distance estimate d « 1.25£’34^kpc. Note 
that the spindown power estimate for an ly = 0.1 heating ef¬ 
ficiency then places the source at 1.5 kpc. Given that Nh is 
below the full Galactic column (at b = —37°), we suspect a 
somewhat higher efficiency, lower power and closer distance. 
The USNO Bl.O proper motion above corresponds to v± = 
100±35(ifepckms“^. MSP binaries te nd to be slower than the 
isolated young pulsar population, with lGonzalez et aTI (1201 Ih 
giving {v±) — 90km s“^. This also supports a relatively 
close distance. 
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The hard power-law X-ray spectrum and strong orbital 
modulation indicate a non-thermal source. The light curve 
with an X-ray maximum flanked by two peaks and the 
observed fx/fj ~ 0.005 are quite typica l of BW/RB 
MSP showing evidence for intrabinary shocks (iRoberts et al.l 
120151) . These double peaks are generally centered around pul- 
sar eclipse at super ior conjunction as for PSR 119593-2048 
(iHuang et al.l |2012h or PSR J2239-0533 (iRomani & Shawl 
120111) . but have al so been seen at inf erior conjunction e.g. 
PSR J2129-0429 (IRoberts et al.llMl) . The local minimum 
between these peaks may be a product of either companion 
eclipse or of beaming along intrabinary shock surfaces. The 
low X-ray lu minosity also s uggests a binary in the pulsar- 
powered state (iLinares 112014l) . 

We close by discussing the most puzzling aspect of 
this system, the as ymmetric double peaked optical light 
curye. Intriguingly, [Li, Halpern & ThorstensenI (1201 4ft haye 
recently measured the optical light curye of the redback PSR 
J1628—3205, which has a 5 h hour period and is near-Roche 
lobe filling. This binary also shows asymmetric maxima and 
minima. These authors do not find a clear origin of this asym¬ 
metry, but speculate on the possibilities of swept back intrabi¬ 
nary shocks or magnetically-directed heating of spots on the 
companion surface. Both are pl ausible for J2049, as well. 
One d ifference, howeyer, is that [Li, Halpern & ThorstensenI 
(120141) place the optical maximum (pB = 0.5 with respect to 
the pulsar ascending node. Thus, with this phasing the star 
is brightest when yiewed at quadrature. Accordingly they 
assign the basic double-peak modulation to ellipsoidal yari- 
ations, with the intrabinary shock or magnetic heating sup¬ 
plying the required asymmetry. In our case without a pulsar 
ephemeris, we do not know the true phases and for the mo¬ 
ment arbitrarily setpB = 0 at the optical maximum. Howeyer 
as noted aboye the characteristic double peaked X-ray light 
curye is often along the pulsar-companion axis. For J2049 
this double peak is in phase with the optical maximum. If we 
interpret that maximum as pulsar inferior conjunction (yiew 
of the heated face), then the double X-ray peak is similar to 
that of PSR J2129—0429. Of course, none of these phase re¬ 
lationships are clear until we obtain a kinematic ephemeris to 
reference our photometric epoch, either through a pulsar dis- 
coyery or through optical spectroscopy. We are pursuing both 
paths. 
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